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Abstract 

We compare QCD theoretical predictions for heavy flavoured me- 
^ ' sons fragmentation spectra in e + e~ annihilation with data from CLEO, 

■ BELLE and LEP. We include several effects in our calculation: next- 

to-leading order initial conditions, evolution and coefficient functions. 
Soft-gluon effects are resummed at next-todeadingdog accuracy. A 
matching condition for the crossing of the bottom threshold in evolu- 
tion is also implemented at next-todeading order accuracy. Important 
initial-state electromagnetic radiation effects in the CLEO and BELLE 
data are accounted for. We find that, with reasonably simple choices 
of a non-perturbative correction to the fixed-order initial condition for 
the evolution, the data from CLEO and BELLE can be fitted with re- 
markable accuracy. The fitted fragmentation function, when evolved 
to LEP energies, does not however represent fairly the D* fragmenta- 
tion spectrum measured by ALEPH. Large non-perturbative correc- 
tions to the coefficient functions of the meson spectrum are needed in 



order to reconcile CLEO/BELLE and ALEPH results. 
Non-perturbative parameters extracted from the fits to e + e~ frag- 
mentation data for D/D* and B mesons are tabulated. They can be 
employed in the theoretical predictions for the production of charmed 
and bottomed mesons in hadron-hadron, photon-hadron and photon- 
photon collisions. 
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1 Introduction 



The study of the fragmentation functions of heavy flavoured hadrons is of 
considerable interest in several aspects of QCD and collider physics. On 
one hand, at transverse momenta much larger than the mass of the heavy 
flavour, the heavy-flavour fragmentation functions behave similarly to the 
light-hadron ones, and obey an Altarelli-Parisi evolution equation. Unlike 
the case of light hadrons, however, heavy-flavour fragmentation is very hard, 
and thus probes a region of the evolution equation in the large x regime. Fur- 
thermore, it is dominated by the non-singlet component, at least at moderate 
energies, so that its study is much simpler. The initial condition for evolution 
has a well defined perturbative expansion in powers of a s (m), m being the 
mass of the heavy flavour, and a well defined expansion in terms of Sudakov 
(i.e. large-x) logarithms. Further corrections of non-perturbative origin, sup- 
pressed by powers of the ratio A QCD /m, can be parametrized to attempt to 
give a uniform description of charm and bottom fragmentation functions. 
These parameterizations can then be employed to provide theoretical pre- 
dictions for charmed and bottomed hadron production in hadron-hadron, 
photon-hadron and photon-photon collisions at large transverse momenta. 

Recently, new, high quality data on charmed meson production have 
come from CLEO [1] and BELLE [2]. One thus has the opportunity to 
perform a more accurate fit to the non-perturbative initial conditions, and 
furthermore one can test the evolution of the fragmentation function from 
centre-of-mass energies of 10.6 to 91.2 GeV, using charm data from LEP ex- 
periments. In the present work we will carry out this program. We develop 
a procedure that overcomes various difficulties in the large-x region of the 
fragmentation function. We are thus able to fit the measured fragmentation 
functions using next-to-leading logarithmic (NLL) evolution, next-to-leading 
order (NLO) initial conditions, NLO coefficient functions, NLL Sudakov re- 
summation (both for the initial conditions and for the coefficient functions) 
and a phenomenological non-perturbative component. With a suitable choice 
of this non-perturbative component we can obtain very good fits to CLEO 
and BELLE data for D* and D fragmentation, over the whole x range. All 
the moments of the fragmentation functions are therefore well reproduced. 
This represents an improvement over previous investigations where, while 
obtaining good fits to some low moments (a necessary and sufficient condi- 
tion for predicting heavy-meson production in hadronic collisions) the fit in 
x-space was not completely satisfactory. The same procedure is also applied 
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to B meson spectra measured in Z° decays. The relevant data from the 
ALEPH [3] OPAL [4], SLD [5] and DELPHI [6,7] collaborations are equally 
well described. 

We evolve the D* fragmentation function fitted to CLEO and BELLE 
data up to LEP energies, taking into account the opening of the bottom 
threshold [8]. We find a discrepancy between the QCD prediction and the 
ALEPH data [9] , that can be parametrized as a power correction of the form 
C(N — 1)/E, with C of the order of few hundreds MeV, or of the form 
C(N - l)/E 2 , with (7 si 5 GeV 2 , where E = ^[q 1 is the total centre-of- 
mass energy and N is the moment in Mellin space. Unfortunately, there 
is no way, at the moment, to discriminate between the two possibilities. 
Theoretical arguments based upon renormalons disfavour the presence of 
1/E corrections in the evolution of fragmentation functions. On the other 
hand, these arguments require validation. 

In Section 2 we describe the theoretical ingredients that enter our calcu- 
lation, collecting and summarizing previously available results: the pertur- 
bative initial condition, the evolution, the Sudakov effects and the bottom 
threshold. The novel treatment of the large-iV/large-a; region is also detailed 
here. In Section 3 we describe our treatment of electromagnetic initial-state 
radiation. The implementation of the non-perturbative component of the 
fragmentation function is discussed in Section 4. In Section 5 we perform 
fits to the CLEO and BELLE data. In Section 6 we compare the evolved 
CLEO/BELLE D* fit to the ALEPH data, and discuss in detail the re- 
lated problems. Fits for B meson production are presented in Section 7. 
In Section 8, we perform fits to data under a different perspective, using 
Mellin moments rather than the x-space distributions, and employing a sim- 
pler, one-parameter non-perturbative function, that can be related to A/m 
power corrections. The implications of the new BELLE and CLEO data for 
heavy-flavour hadroproduction are also discussed in this section. Finally, in 
Section 9, we give our conclusions. 

2 Theoretical framework 

We consider the inclusive production of a heavy quark Q of mass m 

e + e- - Zh(q) - Q (p) + X , (1) 
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where q and p are the four-momenta of the intermediate boson and of the 
final quark. We define x as the scaled energy of the final heavy quark, 



In our framework, we neglect corrections suppressed by powers of the heavy- 
quark (meson) mass, and so the above definition may be replaced with the 
usual experimental definition of scaled momentum (i.e. the heavy flavoured 
meson momentum over its maximum value). The inclusive cross section 
for the production of the heavy quark Q can be written as a perturbative 
expansion in a s 



^(x, q 2 , m 2 ) = J2 *M ^ m 2 , p 2 ) , (3) 

n=0 

where E = a/? 2 is the total centre-of-mass energy, p is the renormalization 
scale, and 

- ^ • (4) 

The cross section (3), normalized to the total cross section, is sometimes 
referred to as the heavy-quark fragmentation function in e + e~ annihilation. 

In order not to spoil the convergence of Eq. (3), the coefficients should 
be small enough to justify a perturbative expansion in terms of a s . There 
are, however, at least two interesting regions of the parameter phase space 
where such convergence is undermined: 

1. If q 2 3> m 2 , large logarithms of the form \og(q 2 /m 2 ) appear in the 
differential cross section (3) to all orders in the perturbative expan- 
sion. These logarithms have collinear origin, and the mass m acts as a 
regulator. 

2. In the region of the phase space of multiple soft-gluon emission, i.e. x — > 
1, the differential cross section contains enhanced terms proportional 
tolog n (l-:r)/(l-a;). 

In the following two sections, we collect the relevant formulae for the 
resummation of these large contributions at the next-to-leading log level. 
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2.1 Collinear logarithms 

In the limit where power terms of the ratio m 2 /q 2 can be neglected, the 
differential cross section satisfies the factorization theorem 

d(jp Q / 9 q x x f 



dx 



x,q 2 ,m 2 ) = Y [ —C Pii (z,q 2 ,fx 2 )D i (-,fj 2 ,mA , (5) 
j x z \ z / 



where the subscript P stands for either T for transverse, L for longitudinal, 
A for asymmetric or nothing for the total (i.e. L+T) cross section 1 . In the 
following, we shall always drop the polarization subscript, since we shall 
always refer to total cross sections. 

The Cj coefficients are the MS-subtracted partonic cross sections for pro- 
ducing the massless parton i, and D; L are the MS fragmentation functions for 
parton i to evolve into the heavy quark Q. The factorization scale /i 2 must 
be taken of order q 2 in order to avoid the appearance of large logarithms of 
q 2 1 \j? in the partonic cross section. The explicit expressions for the partonic 
cross sections and for the fragmentation functions at NLO can be found in 
Refs. [10,11], 

The MS fragmentation functions Di obey the Altarelli-Parisi evolution 



equations 2 



dDl 'x,fx 2 ,m 2 ) = Y Pji ("'"s^ 2 )) Djiz^ 2 ,™ 2 ) ■ (6) 



d log jj? 



The Altarelli-Parisi splitting functions Pji have the perturbative expansion 
PjifaoL.tf)) =a s (fi 2 )P^(x) + a 2 s ^ 2 )p£\x) + 0(a 3 s ) , (7) 



where the P^ ] are 3 [12] 

1 + z 2 3 



1 We follow closely the notation of Ref. [10]. 

2 Notice that the splitting functions are transposed with respect to the structure function 
evolution equations. 

3 The + distribution is defined as 



/ dz h{z) [g(z) ]+ ee [ dz [h(z) - h(l)} g(z) . 
Jo Jo 
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2C 



z 



l-z)+ 

l + a-z) 2 



+ — +2(1-2) + 

z 



11 

12 



n t T B 



5(1 -z) 



(8) 



1 

2 



(9) 



T F [z 2 + (1 - 2) 2 ] , 
n f is the number of active flavours and 

C A = 3, C F = 

The NLO splitting functions Pj^ (needed to achieve NLL accuracy) have 
been computed in Refs. [13-17]. They are too lengthy to be replicated here. 

The initial conditions for the MS fragmentation functions were first ob- 
tained at the NLO level in Ref. [11]. They are given by 

D Q (x, nl,m 2 ) = 5{l-x)+a s {fil)d ( Q ) (x,fil,m 2 ) + 0{a 2 s ) , 
D g (x,nl,m 2 ) = a s (ixl)df\x,iil,m 2 ) + 0(a 2 s ) , 

(all the other components being of order a 2 ), where 

"1 + x 2 



(10) 
(11) 



d%\x,fAQ,m 2 ) = C ¥ 



l-x 



i4 



log^-21og(l-x)-l 



df\x^lm 2 ) = T F [x 2 + (1 - x) 2 ] log 



t4 

m 2 



(12) 



(13) 



In order to compute the NLL resummed fragmentation function, one takes 
the initial conditions at a scale /x — m, evolves them up to \i ~ E (these 
choices for /j,q and /i prevent the appearance of large logarithms that would 
spoil the NLL accuracy), and then applies Eq. (5), using the NLO expression 
for the partonic cross sections Cj given in Eqs. (2.15) of Ref. [10] 4 



C q (z,q 2 ,fi 2 ) = 



5(1 — z) + a s a q 



^ q (q 2 ) , 



C g (z,q 2 ^ 2 ) = 5,0?) 



% 2 



^o, 9 (q 2 ) , 



(14) 
(15) 



4 In this work, we complement the Born clectroweak cross section with a threshold 
factor for the heavy quarks (and antiquarks) 



0-0,9 (<? 2 ) <W<7 2 ) 1 + 



2m 



4m^ 



for q = c,b. Its numerical impact is, however, negligible at the energies considered here. 
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where, to make contact with the notations of Refs. [10, 11], we have defined 

>£) = C T c tlg [z, >£) . (16) 



•Si'- 



The procedure outlined above guarantees that all leading and next-to-leading 
logarithmic terms of quasi-collinear origin (terms of the form (a s log(g 2 /m 2 ))™ 
and a s (a. s log(g 2 /m 2 ))™ respectively) are correctly resummed in the final cross 
section. 

When dealing with the type of convolution appearing in Eqs. (5) and (6), 
it is customary to introduce the Mellin transform 



f(N) = f 
Jo 



dxx N - 1 f(x) . 



(17) 



We adopt the convention that, when N appears instead of x as the argument 
of a function, we are actually referring to the Mellin transform of the function. 
The Mellin transform of the factorization formula (5) is given by 



a Q (N,q 2 ,m 2 ) = £ Q(N, q 2 , /i 2 ) D t (N^ 2 ,m 2 ) , 



(18) 



where 



dxx N ~ 1 — (x,q 2 ,m 2 ) , 



(19) 



and the Mellin transform of the Altarelli-Parisi evolution equation (6) at 
NLO is 



dDj(N, fi 2 ,m 2 ) 
d log fx 2 



(20) 



2.2 Soft logarithms 

Both and d^ contain terms associated to the emission of a soft (and 
collinear) gluon, These terms give rise to a large-iV growth of the corre- 
sponding Mellin transforms 



a^(N,q 2 ,^ 2 ) = C F 



In 2 N + - + 2 lE - 2 In \ In N + a q + O (1/N) 



/// 



(21) 
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m 



-2 ln z iV + 2 ( In — - 2 lE + 1 ) In iV + 8 Q 



where r ) E = 0.5772 ... is the Euler constant and 
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5 2 9 

6 n ~2 +1e+ 2 1e + 



7T 



S Q = 2- — + 2 1e -2 1e 



m 



- - 2 lE In — . 



0(1/N) 
(22) 



(23) 
(24) 



In Ref. [18], the all-order resummation of the large (IniV) contributions has 
been performed to next-to-leading log accuracy, that is, all logarithms of 
the form a™ ln n+1 N (leading logarithms) and a™ ln n N (next-to-leading log- 
arithms) have been correctly resummed. In the following we summarize the 
results of Ref. [18]. 

The Sudakov resummation factor for the e + e~ coefficient function can be 
written as 

A?(JV >(Z V) = exp [lnNgU(\) +g^(\)} , 



where 



S (2) (A) = 



[A + (1 - A) ln(l - A)] , 



2nb 3 



[2A + 21n(l-A) + ln 2 (l-A)j 



(25) 
(26) 

(27) 



+ 



(£« _ 2A^ lE ) 
2tt^ 



ln(l - A) 



and where 
b 



nb 

llC A -4T F n t 



[A + ln(l-A)] 



fi 2 



nb 

nC\ - WC A T F n t - QC F T F n c 



A 4 



12tt -i 24tt 2 

are the first two coefficients of the QCD /3-function, and 



B(1) = A "' 



(28) 

(29) 
(30) 



The variable A is defined by 



\ = b a s (/i 2 ) In TV. 



(31) 



The number of quark flavours in 60 an d b\ is set to the number of active 
flavours at the scale //, i.e. typically four for charm production below or near 
the bottom threshold, and five for charm or bottom production above the 
bottom threshold. 

The presence in the resummed expressions of log(l — A) gives rise to a 
cut singularity starting at the branch point 



This singularity is related to the divergent behaviour of the running coupling 
a s (n 2 ) near the Landau pole at n ~ A QCD , and signals the onset of non- 
perturbative phenomena at very large values of iV or, equivalently, when x 
is very close to its threshold value 1. This translates into an unphysical 
behaviour of the resummed perturbative result in this region. In Section 2.3 
we describe how we have dealt with this issue. 

Similarly to what has been done for the quark coefficient function, in 
Ref. [18] the Sudakov-resummed expression for the initial condition of the 
fragmentation function has also been derived, yielding a result similar to 
Eq. (25). To NLL accuracy we have 




(32) 




(N, & m 2 ) = exp In N gg> (A ) + (A ) 



(33) 



with 



C(Ao) 



AO) 



[2A + (1 - 2A ) ln(l - 2A )] , 



(34) 



27r6 A 




AWb_ 



[4A + 21n(l - 2A ) + ln 2 (l - 2A )] 




and 



= -a 



A = 6 oc s (lJ>l) In N . 



(36) 
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The number of quark flavours in b and bi for the Sudakov resummation 
factor of the initial condition is set to the number of light flavours at the 
scale fio ~ to, i.e. three for charm and four for bottom. Note that, for ease of 
notation, both in Eq. (25) and (33) the renormalization and the factorization 
scales have been taken equal. The full expressions for the Sudakov factors 
can be found in [18]. 

Analogously to Eq. (25), the Sudakov-resummed part Af ni of the heavy- 
quark initial condition also has cut singularities in the complex variable N. 
In the heavy-quark case the singularities start at the branch-point 

' 5 "" (37) 



7V£i = exp 



Ar 



i.e. at A = 1/2 in Eqs. (34) and (35). Again, we defer to Section 2.3 the 
discussion of this problem. 

For later convenience, we introduce here the expansions up to order a s of 
A^ and Af nU defined in Eqs. (25) and (33) 

A*(A,gV) = l + a s ^)[A s q (N, q 2 ,^)] as + 0{a 2 s ) , 

Af ni (7V,^,TO 2 ) = 1 + aM) [Al(N,^,m 2 )] as +0(al) 

where 



(38) 
(39) 



[A s q (N,q*,n 2 )] as =C F 
[Ai i (N^ 2 ,m 2 )] as = C F 



ln 2 A + \- + 2 lE 



fi 2 



In A 



-2 In 2 A + 2 ( In — 



TO 

t4 



2 lE + l) In A 



(40) 



(41) 



Note that they differ from the exact coefficient function of Eq. (21) and from 
the initial condition for the fragmentation function of Eq. (22) only by terms 
finite in the large- A limit. 

In order to merge the NLL-resummed and the NLO expressions without 
double-counting O (a s ) logarithmic terms, we define the Sudakov-resummed 
expressions for the coefficient function and initial condition in the so-called 
'log-R matching scheme' as 

q cs (A^V) = A*(AT,gV) 

xexp{a s (/i 2 ) [a«(iy,g 2 ,/i 2 )- [Af (iV, g 2 , /x 2 )] J } a 0j(? (g 2 ),(42) 

(N,fxl,m 2 ) = Af ni (A,/x 2 ,TO 2 ) 

exp{a s (/i 2 ) [4\N,&m 2 )- [Af ni (A, M 2 , to 2 )]J } . (43) 



D TCS 

Q 



X 
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This matching prescription differs from the one employed in Ref. [18] (see 
Eqs. (36) and (76) there). However, since the exponents in the exponentials 
in Eqs. (42) and (43) are small (i.e. do not contain large logarithms and are 
of order a s ), the exponentials can also be expanded without loss or gain of 
accuracy, giving rise to different - but equivalent - matching prescriptions, 
among which that of Ref. [18]. 

2.3 The large- AT region 

We have previously remarked how the soft-gluon resummation factors and 
Af ni contain singularities at large N which signal the eventual failure of per- 
turbation theory and hence the onset of non-perturbative phenomena. The 
matching of perturbative results with non-perturbative physics is a delicate 
problem, which rests, first of all, on a proper definition of the perturbative 
series. 

One way to address this problem is to work in the framework of infra- 
red renormalons. In Ref. [19], the perturbative series is first improved by 
addition of all sub leading logarithms a" ln fc N, with k < n + 1, in the so- 
called large-/5o approximation. The asymptotically divergent series is subse- 
quently regulated either by truncation at the smallest term or with a Cauchy 
principal-value prescription of its Borel antitransform. This also implicitly 
defines non-perturbative terms which can be cast in the form of power cor- 
rections, hence allowing to relate charm and bottom hadronization. This 
procedure makes maximal use of the insight that can be gleaned from per- 
turbative QCD. However, it will be shown in Section 6 how charm frag- 
mentation data at T(AS) (10.6 GeV) and Z° (91.2 GeV) energies cannot be 
described simultaneously within perturbation theory. Without a more spe- 
cific understanding of the origin of this problem, it would appear premature 
to even attempt to relate rigorously the charm and bottom non-perturbative 
fragmentation functions. 

In the present work, we do not attempt a rigorous formulation of the 
perturbative/non-perturbative matching problem and of the ensuing descrip- 
tion of the non-perturbative terms. We instead simply look for a formulation 
of the resummation prescription that 

(i) is consistent with all known perturbative results, 

(ii) yields physically acceptable results, 
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(iii) does not introduce power corrections larger than generally expected 
for the processes in question, i.e. NA/m for the initial condition [19- 
22] and NA 2 /q 2 for the coefficient functions [23], where A is a typical 
hadronic scale of a few hundreds MeV. 

In detail, as far as the coefficient function is concerned, we make the 
following replacement in Eq. (31) (and hence (25) and (40)) 



1 + fN/Nf 



N -> N j AT /ATT > ( 44 ) 



where N£ is given in Eq. (32) and / is a parameter not smaller than one, 
but of order one. For iV <C N^, the replacement (44) amounts to a tower 
of power corrections to N, starting with f(N — 1)/Ng p» f(N — l)A 2 CD //i 2 , 
consistently with items (i) and (iii) listed above. Furthermore, for large N, 
the replacement (44) becomes N — > N^/f. Thus, with this replacement, the 
functions g^ 1 / 2 ) of Eqs. (26) and (27) have no singularities in the half plane 
Re(iV) > 0, so that item (ii) above is also fulfilled. 

For the initial condition, we apply the same prescription of Eq. (44) 
in Eq. (36) (and hence (33) and (41)), replacing N£ with ATA (defined in 
Eq. (37)) ' 

1 + f IN L - 

N N J 1 1111 (45) 

In this case, the replacement amounts to a tower of power corrections starting 
with f(N - l)/iVA w f(N - l)A QCD //i , and for large N the branch cut in 
Eqs. (34) and (35) is never reached. 

The Landau singularity is regulated if / > 1. For / below 1 the effect of 
the Landau pole should be visible. We plot in Fig. 1 the results of varying the 
/ parameter (/ = 0, i.e. no regulator, 0.5, 1 and 1.5), together with the pure 
perturbative result, without Sudakov resummation. First of all, we notice 
the tiny cusp due to the Landau singularities, located around ps 7.2, 
consistently with Eq. (37). The moments become negative (and therefore 
unphysical) after the cusp. With increasing /, the cusp is displaced to larger 
values of N, until it disappears for / = 1. The fixed order cross section also 
changes sign at A^ ps 11, larger than N^. This is consistent with the large 
N behaviour of Sq shown in Eq. (22), such that Dq(N) becomes negative 
when 

N ~ exp \\ n 2^ • ( 46 ) 
V C F a s {m 2 ) 
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Figure 1: Moments of the perturbative fragmentation function for charm 
production at = 10.6 GeV with and without soft-gluon resummation 
with different values of /, Eqs. (44) and (45), in the regularization of the 
Landau pole singularities. 



However, for small enough a s , this value should be parametrically smaller 
than iV"A. In the case of charm production, this does not happen, remind- 
ing us that we are at the limit of validity of perturbation theory. A more 
consistent behaviour is observed in the bottom case, Fig. 2. In this figure, 
the pattern of improvement towards the large- N region, when going from 
the purely-perturbative initial condition without soft gluon resummation, 
to the inclusion of soft-gluon effects, and then to the addition of the non- 
perturbative regularization, is clearly visible. 

We conclude our discussion with the following remarks. We have found 
that, while for bottom production the NLO result, the inclusion of Sudakov 
effects and the regularization of the Landau singularities follow numerically 
the correct pattern of improvements, in the case of charm production the 
inclusion of Sudakov effects induces a worse large- N behaviour of the cross 
section, signaling the imperfect applicability of perturbation theory in this 
case. Nevertheless, in both cases we have shown that we can obtain a sensible 
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physical result with formulae that are consistent with all known results in 
perturbative QCD, and modest power suppressed effects according to the 
item (iii). We will thus apply our procedure to the fits of charm and bottom 
data. We will use the value / = 1.25 in our fits, since we found that good 
fits can be obtained with this choice. We remark that the parameters of the 
non-perturbative fragmentation function that we obtain in our fits do depend 
upon the choice of /, to an extent that can be inferred from Figs. 1 and 2. It 
is also clear that, for moments around iV pa 5, the effect of the inclusion of a 
regulated Sudakov is modest. We will discuss in Section 8 the implications 
of this fact for hadronic cross sections. 

2.4 The bottom threshold 

In analogy with parton distribution functions, also parton fragmentation 
functions obey matching conditions when crossing heavy-flavour thresholds. 
In Ref. [8], we have computed these matching conditions at next-to- leading 
order in the strong coupling constant a s in the MS scheme. We obtain, 
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neglecting 0(a 2 s ) corrections 



Dq/q( X ' A*thr) m thr) 



*dy 

y 

a s 



D g {x/y,fi 2 hT ,m 2 hr ) 
l + (l-y) 2 



y 



log 



A'thr 



m 



- 1 - 2 logy 



{x, /ithr) m thr) _ ^i/t' ( X > A*thr) m thr) 



2 2 \ 

X, /i t h r , m thrJ 



,.2 



i/i 



i/i 



thr 

3tt ' ta 
for ? 



mr 



(47) 
(48) 



9i,---5n L , (49) 



where n L = n f — 1 is the number of light flavours. Since, in the present 
paper, we are interested in the evolution of charm fragmentation function 
from lower scales, of the order of the charm mass, to higher scales, these 
matching conditions should be used for consistency when crossing the bottom 
threshold. 

In this framework, at low energies (i.e. not much above the charm mass), 
the charm is treated as a heavy quark, in order to provide a perturbative 
expression for its fragmentation function. Near the bottom threshold, the 
bottom is treated as heavy, while all other quarks (including charm) are 
considered light. 



2.5 Simplified evolution scheme 

For the phenomenological analysis performed in the present work, we have 
numerically solved the full set of evolution equations. It turns out, however, 
that, for the case of charm production at T(4S) energies, the contribution 
coming from gluon-splitting processes is fully negligible. Our results, in this 
case, can thus be obtained in a simplified framework, where only the P qq 
splitting function is kept. The Mellin transform of P qq can be performed 
analytically, and one can work with n L flavours, since the annihilation energy 
is of the order of the bottom mass. The evolution equation has the simple 
solution 



E(N,[i ,fi ) = exp <^ log 



a s (n 2 

a s (/Xo) - a s (n 2 ] 
4ir 2 b 



aM) Pi 0) (iV) 



l 27T&0 



2tt6i 



■ (50) 
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Our final formula for the cross section, neglecting singlet contributions, is 
then, using Eqs. (42), (43) and (50) 

a Q (N, q\ m 2 ) = C^N, q\ /i 2 ) E(N, /i 2 , /i 2 ) D™(N, /i 2 , m 2 ) . (51) 

The Mellin transforms for alp, Sp and PjjP are given in formulae (A. 12), 
(A.13) and (A.20) of Ref. [II] 5 . 



3 Electromagnetic initial-state radiation 

Electromagnetic initial-state radiation (ISR) can significantly affect the single- 
inclusive distribution of charmed mesons. At CLEO and BELLE energies, 
the hadronic cross section decreases as the inverse of the squared mass of the 
hadronic system. Initial-state photon radiation is suppressed by a factor of 
a em , enhanced by a logs/m 2 (m e being the electron mass), and, depending 
upon how much energy is radiated away, due to the lower hadronic mass, it 
is enhanced by a larger hadronic cross section. The shape of the fragmenta- 
tion function is also affected, since a consistent fraction of events takes place 
at lower hadronic invariant mass. CLEO and BELLE do not correct their 
data for initial-state radiation, so, in order to perform a meaningful fit to the 
fragmentation function, we have to take it into account. 

We correct the initial distributions of measured inclusive cross sections bin 
by bin, i.e. we find, by an iterative procedure, a new distribution that repro- 
duces the measured one after ISR has been implemented. More specifically: 
calling Xi, % — 1, . . . n, the centre of the bins of the experimental distribution, 
we find a distribution D c {x) (where the suffix c stands for "corrected" ) that 
is continuous, vanishes at x = and x — 1, and is linear in all intervals 
(0, xi), (xi, X2), ■ ■ ■ , (x n -i, x n ), (x n , 1), such that, when ISR corrections are 
applied, we reproduce the measured distribution. 

We model ISR in the following way. We assume for the radiated electro- 



5 We point out that there is an obvious misprint in Ref. [11], where formula (A. 7) should 
be replaced by 

j»+ib g r(s) 
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magnetic energy the distribution [24-27] 



where 

s 



/3 = 2 ae,n 



log' 2 



— = 5/3(1 -z)>^ ->-(! + z), (52) 



5=1 + / + — (y-2 , z = — , (53) 



7T 

s = g 2 being the squared centre-of-mass (CM) energy, and Shad the square 
of the invariant hadronic mass. The kinematic distribution of the hadronic 
system is assumed to be as if only a single photon, collinear to either the 
electron or the positron, was radiated. This assumption neglects double 
radiation, which gives effects of the order (3 2 , and the transverse momentum 
of the radiated photon, which is typically much smaller than the available 
energy. We use the Born cross section for the heavy-quark production in the 
hadronic reference frame. The value of x in the laboratory frame is obtained 
by a Lorentz boost. In summary 

D ( x i) = f 4m 2 dz j d y d cos 6 a ^ ^ a °^ C o S C fl S ^ ^ Dc ^ 6 ( Xi ~ x ( z >y> 9 )) ■ 

(54) 

where x(z, y, 8) is the momentum fraction of the heavy flavoured hadron (of 
mass m h ) in the e + e~ CM frame. It is obtained as follows. We define the 
momentum components of the hadron in the hadronic CM system 

Ph = ^sz-4m 2 h , p° h = ^p 2 h + m 2 h , pl=p h cos6, (55) 

so that y is its momentum fraction. Then we boost it to the e + e~ CM frame. 
Under our assumptions (that all the electromagnetic energy is collinear either 
to the electron or to the positron, and that double radiation and the photon 
transverse momentum are negligible) the boost velocity is purely longitudi- 
nal, and is given by v — (1 — z)/(l + z). The hadron momentum in the e + e~ 
CM frame is then 

o_Ph + v p\ 



and 



V 



P° = ^=^, P=\hl-<, (56) 



X (z,y,e) = -=1== T . (57) 
V s / 4 - m h 
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We use for o"o the exact Born cross section in the massless limit. Since 
asymmetries cancel in Eq. (54), we always assume the angular dependence 



1 da (zs, cos ( 
<Jq(s) dcos6 



a (zs) 3 , 2 
\ / - (1 + cos 2 60 . 



(58) 



We do, however, supply the threshold factor to the total cross section. Thus, 
near the T(45) we have 



cr {zs) 
o-o(s) 



, i i 2m 2 

6{zs - Am 2 h ) ^ ~ 



Am 2 

sz 



Z 



4m 2 



(59) 



We checked that the effect of finite mass corrections to the angular dis- 
tribution, the use of mt instead of the quark mass m in the threshold factor, 
as well as the scaling violations in the fragmentation function due to the 
reduced (i.e. s — > sz) CM hadronic energy, have a negligible impact on our 
results. 
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Figure 3: Left: the effect of the ISR correction on BELLE data for D* + — ► 
D°ii + . Right: the same data in moment space, shown together with the 
ALEPH ones. 



The effect of the ISR correction is displayed in Fig. 3, where BELLE data 
for D* + — ► D°7T + are displayed before and after the ISR correction has been 
applied, both in x and moment space. The ALEPH data are also shown in 
moment space. We see that the corrected spectrum for BELLE is harder 
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and lower in normalization than the uncorrected one. This is to be expected, 
since ISR lowers the available hadronic energy, thus softening the spectrum 
and increasing the cross section at the same time. We can also see that the 
effect for BELLE is not large, but nonetheless not negligible. It is instead 
much less prominent, up to the point of being negligible, for the ALEPH 
data taken on the Z° peak, as expected. 

4 Non-perturbative fragmentation function 

In the heavy-quark fragmentation-function formalism, the largest non-perturbative 
effects come from the initial condition, since one expects power corrections 
of the form A/m. We assume that all these effects can be described by a 
non-perturbative fragmentation function D NP , that takes into account all 
low-energy effects, including the process of the heavy quark turning into a 
heavy-flavoured hadron, that has to be convoluted with the perturbative 
cross section. Thus, the Mellin transform of the full resummed cross section, 
including non-perturbative corrections, is 

a H (N,q 2 ) = a Q (N, q 2 , m 2 )_D NP (iV) . (60) 

We have attempted to fit CLEO and BELLE D* data using several forms for 
Dnp- We found that the best fits are obtained with the two-component form 

D NP (x) = Norm, x -±- [5(1 - x) + cN^(l - x) a x b ] , (61) 

with 

N a , b = [ (l-x) a x b . (62) 
Jo 

This form is a superposition of a maximally hard component (i.e. the delta 
function) and the form proposed in Ref. [28]. It can be given a simple phe- 
nomenological interpretation, the hard term corresponding in some sense to 
the direct exclusive production of the D*, and the Colangelo-Nason form 
accounting for D*'s produced in the decay chain of higher resonances. 

Following the approach of Ref. [29], we assume that the D meson non- 
perturbative fragmentation function is the sum of a direct component, which 
is isospin invariant, plus the component arising from the D* decay. The decay 
D* — > Dn is very close to threshold, so that the D has the same velocity 
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of the D*, and their momenta are thus proportional to their masses. Under 
these circumstances, the component of the D fragmentation function arising 
from D* — > Dtt decays is given by 

B(D* -> Dn) D°{x) , (63) 

where we have defined 

tt (x) = D NP x—— — — 0(1- x—— , (64) 
\ m D J m D \ m D J 

and B(D* — > Dn) is the branching ratio of D* — ► Z?7r. Observe that Z)^ has 
been defined so as to have the same normalization as D® P . In N space we 
obtain immediately 



D?(N) = D&(N) 



N-l 



m D * 



(65) 



For the D* — > D7 decay, in the £)* frame, the £) has non-negligible 
velocity, but it is non-relativistic. We call 9 its decay angle with respect to 
the D* direction, and pa its momentum 

2 2 

m D* ~ m D / aa \ 

Pd = — • 66 

2m D * 

We call (5 the D* velocity and 7 = — (3 2 . Thus, the longitudinal 

component of the D momentum in the laboratory frame is given by a Lorentz 
boost 

7 (p D cos 6 + j3m D ) , (67) 

where we have neglected terms of order p 2 D . Thus the component of the D 
fragmentation function coming from D* — > D7 decay is given by 

B(D* -> D 7 ) ^(x) , (68) 

with 

b ?(x) =[ d y d -^l DSM S ( ■<<*•>«»» + fro) _ \ _ {m) 

J 2 V Pmax / 

where j9 max is the maximum D momentum in the laboratory. Since we always 
consider the ultra relativistic limit, we have 

V = ^L, /J-l, (70) 

Pmax 
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so that we obtain 



y(p D cos6 + m D ) 



x 



(71) 



The double integral cannot be performed in closed form. However, it is easy 
to obtain the moments 



D°(N) = D&(N) 



dcos8 



Pd cos 6 + rriD 
m D * 

N 



N-l 



m D * (m D + p D ) - ( m D - Pd) 
2Pd 



(72) 



We thus describe D + /° production as the sum of a primary (i.e. not coming 
from D* decays) component, plus the contributions coming from D* decays 



D® p (x 



D°t' P (x) + B(D*+ -> D + 7T°)D 



\ hD+, 



NP K-^J -r ' ^ " \X) 

+B(D* + -> D + 1 )D^ + D{x) , (73) 
NPK x) = DZ' P (x) + lB(D*+^D°n+) + B(D* ^D°n )]Df(x) 

+B(D*°^D° 1 )D°°(x) . (74) 



D D ° 



We took the value of masses and branching ratios from Ref. [30]. For ref- 
erence, we report in Table 1 the values we used for the masses and for the 
decay rates of the charmed mesons. 



D (mass in GeV) 


branching ratios 


£>*°(2006.7±0.4) -> DV 
— > D°7 


0.619 ±0.029 
0.381 ±0.029 


£>*+(2010.0±0.4) -> D°tt+ 
-> D + 1 


0.677 ±0.005 
0.307 ±0.005 
0.016 ±0.004 


J D°(1864.5 ± 0.4) -> K-71+ 
L>+(1869.3 ± 0.4) -> K-71+71+ 


0.0381 ±0.0009 
0.092 ±0.006 



Table 1: Charm hadron masses and branching ratios. 
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5 D mesons data fits near the T(4S) 



Several parameters enter our calculations. First of all, at all matching points, 
there are scale choices that could be varied, to yield a perturbative uncer- 
tainty in our result. Those are the initial evolution scale /i , the matching 
scale for the crossing of the b threshold // t h r , and the final evolution scale \x. 
In the present work we fix 

(J, = m, fi= \fc? , /ithr = m th r = m b . (75) 

These scales could, in principle, be varied by a factor of order two around 
the values listed above, yielding a sensibly different result. However, in gen- 
eral, the scale variation will simply result in different values for the fitted 
parameters of the non-perturbative form. When computing cross sections 
for different processes, one should then use the parametrization appropriate 
for the scale choice that has been made in the fit, hence compensating for the 
change. In the present work we will not pursue this issue further, since our 
aim is simply to show that a fit within QCD is possible. A similar remark 
applies to the value of A QCD and the quark masses, that we will fix at 

Aq^ d = 0.226 GeV , m c = 1.5 GeV , m b = 4.75 GeV . (76) 

The CLEO and BELLE data are given as absolute cross sections. Since 
we correct the data for ISR effects, we should normalized our data to the 
e + e~ charm cross section corrected for ISR effects. We thus use the value of 
R(e + e~) reported in Ref. [31], defining 

where (s) = 86.86 nb/s is the Born cross section for e + e~ — > 

3.56 is the value of R measured by CLEO, 0.4 is the charm fraction, and the 

factor of 2 allows for the counting of both charge conjugate states. 

We have fitted all D* + and D*° data with the same set of parameters, 
except for the normalization, which is kept independent for each data set. 
This procedure is justified, since the errors in the data do not include overall 
errors that do not affect the shape of the fragmentation function. We have 
limited ourselves to the fit range 0.2 < x < 1 for CLEO and 0.08 < x < 1 
for BELLE. In the case of BELLE data, we use only the continuum sample 
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for x < 0.5, and for x > 0.5 we combine the continuum and the on-resonance 
sample in the following way 

yjs 2 c + y r /s 2 r 1 d c /s 2 c + d r /s 2 r 

y I/** + 1/* ' + W ' V^ + iAv 2 ' 1 j 

where y c / r , s c / r and d c / r are the central value, the statistical error and the 
systematic error of the continuum/on-resonance data, and y, s and d are our 
combined central value, statistical error and systematic error. For all data 
sets we combine the statistical and systematic errors in quadrature. 



Eq. (61): a = 1.8 ±0.2, b = 11.3 ±0.6, c = 2.46 ± 0.07, total x 2 = 139 


Set 


(C) D*+ 


(B) D*+ -> D° 


(B) D*+ -> D+ 


(C) D*° 


(B) D*° 


Norm. 


0.238 


0.253 


0.227 


0.225 


0.211 


X 2 /pts 


33/16 


63/46 


13/46 


13/16 


17/46 



Table 2: Results of the fit to D* CLEO (C) and BELLE (B) data. The last 
line reports the x 2 over the number of fitted points for each data set. 



The result of the fit is reported i 
we show the data and the fitted curv 
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Figure 4: Fit to 



. Table 2 and in Figs. 4, 5, 6, 7 and 8 
, both in x and moment space. 
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x N 

Figure 5: Fit to BELLE D* + -> D° data. 



A considerable part of D's are produced indirectly through D* decays. 
Here we assume that both D*'s and the D's that are not the product of D* 
decay are produced with a charge-independent rate. Under this assumption, 
the fraction of "direct" (meaning not arising from D*) and "indirect" (from 
D*) D mesons are in relative proportion of 0.473 to 0.527. These numbers 
can be extracted from the total production cross section of charmed me- 
sons reported in Ref. [1], and from Table IX of Ref. [2]. We then use the 
parametrization of Table 2 for the D* production, the branching ratios for 
D* — > D decays given in Table 1 and a description of the decay as detailed 
in Sec. 4. 

We parametrize the direct D component with the same form used for the 
D*, and fit it to the D + production data, where a larger fraction of direct D 
is expected. We then use the fitted direct D parametrization to describe the 
direct part of the D° production data. In all cases, the overall normalization 
is chosen for a best fit to each data set, in order to be insensitive to overall 
normalization differences. 

The result of the fit for the D + /° mesons is reported in Table 3. In Figs. 9, 
10, 11 and 12 we show the data and the fitted curve, both in x and moment 
space. 
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Eq. (61): a = 1.1 ±0.1, b = 7.6 ± 0.6, c = 4.6 ± 0.2 




total x 2 = 50 


total x 2 = 109 


Set 


(C) D + 


(B) D+ 


(C) D° 


(B) 


Norm. 


0.263 


0.270 


0.609 


0.598 


X 2 /P ts 


14/16 


36/46 


32/16 


77/46 



Table 3: Results of the fit to D CLEO (C) and BELLE (B) data. The fit 
was performed over the D + data only, that are more sensitive to the direct 
component, and then used to describe D° data. 
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Figure 7: Fit to CLEO D*° data. 
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Figure 8: Fit to BELLE D*° data. 



25 





Figure 10: Fit to BELLE D + data. 
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6 D mesons data fits on the Z° 

In this work we use the data from the ALEPH collaboration [9], which are 
the most precise ones. These data are affected by electromagnetic initial- 
state radiation as well. However, unlike the T(AS) case, the ISR does not 
appreciably distort the spectrum, but it mostly affects the total cross section. 
This is easily understood: any appreciable amount of ISR on the Z° peak 
brings the reaction off resonance, to a vanishing cross section. Thus, the bulk 
of heavy flavour production always takes place on the Z° peak. Conversely, 
at the T(AS) energy, the ISR generates processes with higher cross section 
and a lower hadronic invariant mass. Since the ALEPH data are normalized 
to the total number of hadronic events, the effects of ISR largely cancel in 
the ratio. We shall anyway perform the correction for ISR for these data as 
well. 

In Fig. 13 we display our fit with ALEPH data. We fit the data in the 
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Figure 13: ALEPH D* + data and the result of our non-singlet fit (solid line). 
The dashed line represents the result obtained using full evolution. 



region x G [0.4, 1] using the non-singlet component only, since a subtraction 
of the gluon-splitting contributions was performed by ALEPH. Observe that, 
in this calculation, the bottom-threshold crossing has to be dealt with, ac- 
cording to the discussion of Section 2.4. We also show, for comparison, the 
full evolution result (dashed line), using the same parameters obtained in the 
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non-singlet fit. As we can see, the difference is only visible at small x. The 
result of the fit for the non-perturbative parameters is 



2.4 ±1.2, b =13.9 ±5.7 c = 5.9 ±1.7, 



(79) 



with a x 2 — 4-2 for 13 fitted points. These results are not really consistent 
with those for the T(4S) data in Tab. 2. 

In order to better quantify the discrepancy between Eq. (79) and Tab. 2 
we use the parametrization of CLEO and BELLE data to predict the D* 
fragmentation function at LEP energies. The LEP prediction, using the 
parametrization of Table 2, is reported in Fig. 14 together with ALEPH data. 
We observe that the fitted normalization is very close to the CLEO D* + nor- 
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Figure 14: ALEPH D* + data, compared to the QCD prediction. 



malization. We find a x 2 — 60.1 (for 13 fitted points) for this parametrization. 
Thus, the description is not satisfactory, especially in the large-x (large-iV) 
region. 

In Fig. 15 we show the ratio of the moments of ALEPH D* + data over our 
prediction. We observe that the N dependence of the ratio is well described 
by the functional form 

1 ±0.044 (N — 1) ' ^ 
where, since the first moment of the non-singlet distribution should be exactly 
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given by the theory (because of charge conservation), we normalize to one 
the extrapolation of the data to N = 1. 

We can only speculate about the possible origin of the discrepancy and 
the form of the coefficient of (N — 1) in Eq. (80). Assuming that we are 
dealing with a non-perturbative correction to the coefficient function of the 
form 

l + (81) 
this would lead to the extra factor 

, , C(N-l) 

1 H W 2 — 

(82) 



1 _l_ C(N-1) ' 

(where Mz and My are the Z° and T(45) mass) to be applied to our pre- 
diction for the ALEPH data. For C = 5 GeV 2 we reproduce the behaviour 
of Eq. (80). In Ref. [23], on the basis of a calculation of infra-red renormalon 
effects, a 1/q 2 power correction is found, with an iV dependence marginally 
compatible with (81). No 1/E correction is found. Ref. [32] also predicts a 



30 



leading 1/E 2 power correction. However, the C ~ 5 GeV 2 coefficient would 
appear to be somewhat too large 6 . Alternatively, if we admitted the existence 
of corrections to the coefficient functions of the form 

l + ^>. (83, 

then we would find C ~ 0.52 GeV, a much more acceptable value. We 
observe that a form 

n C\ N ~ 2 1 C(N-2) 

1 + -J ^1 + ^-^ (84) 

was required in Ref. [10] to fit light-hadron fragmentation data. 

Demonstrating the absence (or the existence) of 1/E corrections in frag- 
mentation functions would be a very interesting result, since it would help 
to validate or disprove renormalon-based predictions. Unfortunately, the low 
precision of the available data does not allow, at the moment, to resolve this 
issue. 

We would like to remark that the discrepancy between the CLEO /BELLE 
and ALEPH data exclusively depends upon the evolution between the T(4S) 
and Z° energies. The method we used to describe the CLEO/BELLE data 
(i.e. the perturbative calculation of the fragmentation function, the Sudakov 
effects in the initial conditions and the parametrization of the non-perturbative 
part) does not affect the conclusions of the present section. In fact, we can 
simply compute the ratio of the moments of the inclusive D* + (ISR cor- 
rected) distribution at CLEO/BELLE and ALEPH, and compare it to the 
theoretical prediction. The result of this comparison (where we have used, 
for simplicity, BELLE data only) is displayed in Fig. 16. The curves are 
given by 

a Q (N,Mlm 2 ) _ a q (N,M%, 2 2 1 + a s (fi 2 r ) / n 

a Q (N,M 2 ,m 2 ) I + c^D/tt ^ iV '^'^ a q (N,M 2 ,^ r ) 

where \xz and /ix are the factorization scales and the evolution factor E is 
given in Eq. (50). Notice that low-scale effects, both at the heavy quark mass 
scale and at the non-perturbative level, cancel completely in this ratio, mak- 
ing its prediction entirely perturbative. For a q , in the NLO results (dashed 

6 If we believe that it is the maximum meson energy, not E, that controls power effects, 
than we would have C w 1 GeV 2 , a more acceptable value. 
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N 

Figure 16: The ratio of ALEPH and BELLE moments for the D* + fragmen- 
tation function, compared to QCD evolution. The solid band is obtained 
with QCD NLO evolution and Sudakov effects in the coefficient functions, 
while the dashed bands is NLO evolution only. The bands are obtained by 
setting nz/r — £,M z /y and varying 1/2 < £ < 2. 



lines), we have used 



a q (N, q 2 , fi' 2 ) = l + a s (fi 2 ) a« (N, q 2 , ^ 2 ) 



(86) 



while for the full result (solid lines) we have included the NLL resummation 
of soft gluon emission in the coefficient functions 



J7) 



a q (N,q 2 ,fi 2 ) = A s q (N,q 2 ,fi 2 ) 

x |l + « s (/i 2 ) [aU(N,q 2 ,n 2 )-[A s q (N,q 2 ,^)] 

The definitions of ctq and are given in Sections 2.1 and 2.2. We have 
set Hz/r — £M z /r with £ = 0.5, 2 to plot our bands. As we can see from the 
figure, the rather large scale uncertainty displayed by the NLO result is much 
reduced when Sudakov effects are included. In both cases, however, the data 
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clearly undershoot the pure QCD prediction, being instead compatible with 
the inclusion of the correction factor (80) (dotted lines). We have also checked 
that our full result is essentially unchanged if, instead of formula (87), we use 
the fully exponentiated formula (42). Furthermore, the change of variable 
given in Eq. (44) to deal with the Landau pole has very little impact on 
our curves. Using the very large value Aqc D = 0.3 GeV would lower the 
theoretical predictions by no more than 11% for iV < 20, very far from 
explaining the observed effect. 

The deconvolution of ISR effects, that hardens the T(45 l ) data, but is 
insignificant on the Z°, widens the discrepancy. However, if we did not 
apply the deconvolution, the effect would still be partially visible. 

Because of the relatively low energy of the data on the T(4S f ), it is le- 
gitimate to wonder whether charm-mass effects could be responsible for the 
discrepancy between LEP and T(45 l ) data. We have not included mass ef- 
fects in the present calculation. However, in Ref. [33], mass effects in charm 
production on the T(4S) where computed at order a|, and found to be small. 
We thus believe that it is unlikely that mass effects could play an important 
role in explaining this discrepancy. 

7 B mesons data fits on the Z° 

The same framework that yields good fits to D meson production data can 
also be used to describe B meson production on the Z°. Accurate data have 
been published by the ALEPH [3], OPAL [4] and SLD [5] Collaborations. 
Preliminary data are available from DELPHI [6,7]. We find that, to describe 
B production, the 5(1 — x) term in Eq. (61) is in fact not needed, i.e. the c 
parameter tends to become very large in the fitting procedure. In this limit 
the form of Eq. (61) becomes a two parameter form, coinciding with that of 
Ref. [28]. In Fig. 17 we show the result of a simultaneous fit to ALEPH and 
SLD data. In Fig. 18 we show the same best-fit curve together with the OPAL 
and DELPHI data. The fit yields a = 24 ±2 , b = 1.5 ±0.2 with a x 2 = 43 for 
SLD and 51 for ALEPH for 21 and 19 data points respectively (for these data, 
bin-to-bin correlations provided by the experimental Collaboration were also 
taken into account). We did not attempt to fit the OPAL data together, since 
it was not clear to us how to handle the asymmetric, correlated systematic 
errors given by OPAL. However, it is clear from Fig. 18 that also this data 
set, as well as the preliminary data from DELPHI, is well described by the 
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Figure 17: Fit of the fragmentation function for B production together with 
ALEPH (upper) and SLD (lower) data. 



fit 7 . 



8 Moment- space fits and power corrections 

The fits presented so far have been performed on the measured x-space dis- 
tributions, and they were aimed at providing an accurate description of all 
the experimental data. This has required a flexible parametrization for the 
non-perturbative fragmentation function, leading to the choice of the three- 



7 Note that, while the ALEPH set refers specifically to B mesons, the SLD, DELPHI 
and OPAL data are for all weakly decaying 6-flavoured hadrons. The two quantities could 
therefore be slightly different, due to the small fractions of B s and B baryons (10% each, 
see Ref. [30]). For an example of a quantitative estimate see Eq. (5.10) of Ref. [7]. 
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Figure 18: The fragmentation function fitted to ALEPH and SLD B data 
shown here together with OPAL and DELPHI data. 



parameter form given in Eq. (61). The data are fitted well in the large-x 
region, so that all moments of the fragmentation function are also well re- 
produced. This is important, since, as noted in Refs. [34,35], heavy-flavour 
production spectra in hadronic collisions are determined by a few Mellin mo- 
ments (usually in the range AT = 2, ... 6) of the non-perturbative fragmen- 
tation function. This property was successively exploited in Refs. [29,36-38] 
for predicting bottom and charm spectra in pp collisions. Inaccuracies in the 
description of the large-x region in e + e~ annihilation could therefore lead 
to large errors in the moments that are relevant to the hadropro duct ion of 
heavy quarks. Conversely, in the framework of heavy-flavour production, an 
accurate fit in x-space is unnecessary, as long as the moments are well fitted 
in the relevant range. For this purpose, it is therefore convenient and suffi- 
cient to use for the non-perturbative fragmentation function one-parameter 
functional forms that are commonly found in the literature [39-41]. In the 
following discussion, we will focus upon these one-parameter forms, and in 
particular on the one of Ref. [39] 

D NP (x) = (a + l)(a + 2)x Q (l-x). (88) 
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Figure 19: Moments of the non-perturbative component Dnp(N) extracted 
from e + e~ D* data, and those of the fitted non-perturbative fragmentation 
function (61) with the parameters of Table 2. 



It is important to stress that the choice of a specific parametrization like 
this one is exclusively a matter of convenience, aimed at easing the transfer 
of the non-perturbative information from e + e~ collisions to other processes. 
One can either choose a different functional form, or even analyze the data 
in terms of non-perturbative moments Dnp(N) only. In Fig. 19, we show the 
moments D-^p(N) extracted from e + e~ data. The points in the figure are 
obtained by taking the experimental values of the moments of the fragmen- 
tation function together with their errors, divided by the pure perturbative 
component of the fragmentation function, computed with our default pa- 
rameters 8 . In the figure we also show the non-perturbative component given 
by the form (61), with the parameters taken from Table 2. Also evident is 
the poor consistency between values obtained from data taken on the T(4S*) 
and on the Z°. This is, of course, the same situation already observed in 

8 The perturbative fragmentation function for ALEPH is computed using the non-singlet 
component only, since gluon-splitting contributions have been subtracted from the pub- 
lished experimental distribution. 
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Figure 20: Fits to D* + data for the parameter a of the parametrization (88). 



Using the Mellin transform of formula (88) 

(q + l)(q + 2) 
Dnp(A ° " (a + iV)(a + iV + l) ' (89) 

we can translate the moments in Fig. 19 into values for a with the appro- 
priately propagated error. The results are displayed in Fig. 20. From the 
figure we see that the one-parameter form (88) does not describe perfectly 
the whole shape, as shown by the non-constancy of a extracted from differ- 
ent moments. However, to a good degree of approximation a single value 
of a can describe all the moments up to N ~ 6 or so. This is enough for 
the purpose of using the fitted function for convoluting a pr distribution in 
hadronic collisions. 
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8.1 Scaling property: from D to B mesons 



Several theoretical arguments in QCD [19-22] predict for the heavy quark 
non-perturbative fragmentation function the behaviour 



where A is a hadronic scale, and m is the mass of the heavy quark. If _D NP (iV) 
depends upon a single parameter, its value can be linked to the ratio A/m. 
For example, in the case of the form (89) the series expansion in powers of 
1/a is given by 



Reinterpreting a — > 2 m/ A we will be able to check the behaviour of the 
leading power correction. Using A ~ 300 MeV, m c ~ 1.5 GeV and m^ ~ 
4.75 GeV, one expects to find old ~ 10 and «b ~ 30 when fitting D/D* 
and B mesons respectively. Whatever the exact values are, it will always be 
possible to test for the predicted scaling law 



To this end, we extract the value of a for B meson production at Z Q 
energy. In Fig. 21 we show the fits to the four available data sets. All 
the data appear consistent with each other. Within fairly large uncertainties 
(resulting from the non-constancy of a through the fits to different moments, 
and the discrepancy between the determination of old at Mx and Mz) we 
can see that the expectations are largely fulfilled, leading to an ag/«o ratio 
of order 1.5 to 3. The tendency for values smaller than mb/m c = 3.17 might 
be a consequence of a number of factors, like the B data being for "weakly 
decaying" _B's, and therefore generally softer than the leading B* and B**, 
or the mass entering the power corrections being closer to the meson mass 
rather than the quark mass. 

It is also worth noting that, given a value for olb — 25, the expectations 
for the value of the ratio are much better fulfilled if we use the old fitted 
at the Z energy {p. D — 9) rather than the one fitted at the T(AS) energy 
(old — 14), as shown in Fig. 22. This result mildly supports the view that 
large non-perturbative corrections may affect the T(45') data. 



D NP (N) = 1 — (N — l)A/m + 0(A 2 /m 2 ) , 



(90) 




(91) 



(92) 
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Figure 21: Fits to weakly decaying 5's data for the parameter a of the 
parametrization of Ref. [39]. 



The results of Figs. 19 and 20 are also summarized in Table 4, together 
with similar results for the one-parameter forms of Ref. [41] and with the 
popular PSSZ form [40] 9 . Results for B mesons are also shown. 

9 The values of e that we find in this case are about one order of magnitude smaller than 
those usually extracted from Monte Carlo simulations. They lead, therefore, to a harder 
non-perturbative fragmentation function and hence to larger rates in hadronic collisions. 
It is moreover worth noting that the use of the PSSZ fragmentation function in a context 
where non-perturbative corrections are expected to scale like 1— (TV— l)A/m is inconsistent. 
In fact, while the e parameter can be interpreted as being of order A 2 /m 2 [40], the series 
expansion of its Mcllin transform can be shown to be 

^ N p(iV) = 1 + 2(log £ -l)(^-l) + 4^(0) W+7g)v - + 

The presence of the log e term in the coefficient of the term linear in y/e does not allow to 
interpret it as a simple A/m power correction. 
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JpQCD 








c ® 10.58 GoV 


0.7359 


0.5749 


0.4601 


0.3778 


0.3167 


0.2698 


0.2331 


c ® 91.2 GeV (NS) 


0.5858 


0.3937 


0.2843 


0.2151 


0.1683 


0.1345 


0.1107 


c @ 91.2 GoV (full) 


0.5954 


0.3988 


0.2860 


0.2158 


0.1686 


0.1353 


0.1108 


6 ® 91.2 GoV 


0.7634 


0.6280 


0.5309 


0.4590 


0.4033 


0.3587 


0.3222 


Experimental data (norm, to one) 


BELLE D* + -> D" (ISR corr.) 


0.6418 ± 0.0042 


0.4399 ± 0.0028 


0.3169 ± 0.0020 


0.2375 ± 0.0015 


0.1838 ± 0.0012 


0.1462 ± 0.0010 


0.1189 ± 0.0009 


ALEPH D'+ (ISR corr.) 


0.4920 ± 0.0152 


0.2803 ± 0.0075 


0.1748 ± 0.0047 


0.1160 ± 0.0033 


0.0806 ± 0.0025 


0.0582 ± 0.0020 


0.0432 ± 0.0016 


ALEPH B 


0.7163 ± 0.0085 


0.5433 ± 0.0097 


0.4269 ± 0.0098 


0.3437 ± 0.0096 


0.2819 ± 0.0094 


0.2345 ± 0.0091 


0.1975 ± 0.0087 
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Tab. 2 and Eq. (61) 


0.868 


0.767 


0.688 


0.626 


0.576 


0.536 


0.503 


ALEPH B 


0.938+^14 


o-ses+o'lils 


"• 8 '' 4 ±0.020 


''• 74 9±0.023 


-f399io.O25 


0.654±' r ; : ^ 4 


0-6131^.029 


SLD B 


o.93il»; ^ 


o.ssotg;"^ 




°- 718 to.S24 


0.661+gg24 


°- 6i °tn24 




KLP a 


CLEO D*+ 


13.28 ± 1.27 


12.76 ± 0.54 


13.21 ± 0.40 


13.70 ± 0.34 


14.23 ± 0.32 


14.81 ± 0.32 


15.43 ± 0.33 


BELLE £>*+ -» D° 


12.64 ± 0.70 


12.49 ± 0.35 


13.16 ± 0.27 


13.82 ± 0.24 


14.51 ± 0.23 


15.24 ± 0.23 


16.01 ± 0.23 


ALEPH D*+ 


9.49 ± 2.03 


9.33 ± 0.93 


9.42 ± 0.68 


9.58 ± 0.59 


9.76 ± 0.56 


9.96 ± 0.57 


10.18 ± 0.59 


ALEPH B 


29.42 ± 5.82 


25.12 ± 3.41 


24.56 ± 2.90 


24.22 ± 2.67 


24.00 ± 2.58 


23.86 ± 2.56 


23.75 ± 2.57 


ALEPH B, m b = 4.5 GeV 


34.32 ± 7.77 


28.26 ± 4.23 


27.34 ± 3.51 


26.73 ± 3.18 


26.34 ± 3.04 


26.07 ± 2.98 


25.87 ± 2.97 


ALEPH B, m b = 5.0 GeV 


25.90 ± 4.59 


22.72 ± 2.84 


22.41 ± 2.46 


22.23 ± 2.30 


22.13 ± 2.25 


22.07 ± 2.24 


22.03 ± 2.26 


SLD B 


25.87 ± 8.66 


22.11 ± 3.80 


21.33 ± 2.75 


20.73 ± 2.25 


20.27 ± 1.98 


19.89 ± 1.82 


19.57 ± 1.72 


BCFY r 


CLEO -D* + 


0.0531 ± 0.0077 


0.0610 ± 0.0036 


0.0615 ± 0.0026 


0.0611 ± 0.0021 


0.0601 ± 0.0019 


0.0587 ± 0.0017 


0.0569 ± 0.0016 


BELLE £>*+ -> D° 


0.0570 ± 0.0046 


0.0628 ± 0.0025 


0.0618 ± 0.0018 


0.0604 ± 0.0014 


0.0585 ± 0.0013 


0.0564 ± 0.0011 


0.0541 ± 0.0011 


ALEPH D'+ 


0.0849 ± 0.0247 


0.0936 ± 0.0125 


0.0972 ± 0.0092 


0.0988 ± 0.0080 


0.0993 ± 0.0074 


0.0990 ± 0.0072 


0.0981 ± 0.0072 


ALEPH D" + , m c = 1.3 GeV 


0.0470 ± 0.0238 


0.0557 ± 0.0102 


0.0594 ± 0.0074 


0.0613 ± 0.0063 


0.0621 ± 0.0059 


0.0623 ± 0.0057 


0.0619 ± 0.0057 


ALEPH D*+, m c = 1.7 GeV 


0.1198 ± 0.0289 


0.1288 ± 0.0146 


0.1323 ± 0.0108 


0.1336 ± 0.0093 


0.1337 ± 0.0086 


0.1329 ± 0.0084 


0.1315 ± 0.0084 


PSSZ e (xlO^) 


BELLE £>* + -» D u 


0.234 ± 0.032 


0.271 ± 0.019 


0.260 ± 0.013 


0.246 ± 0.010 


0.230 ± 0.009 


0.213 ± 0.007 


0.197 ± 0.007 


ALEPH D'+ 


0.473 ± 0.245 


0.548 ± 0.129 


0.574 ± 0.096 


0.580 ± 0.081 


0.575 ± 0.074 


0.563 ± 0.071 


0.547 ± 0.069 


ALEPH B 


0.028 ± 0.014 


0.047 ± 0.016 


0.056 ± 0.016 


0.062 ± 0.017 


0.068 ± 0.018 


0.073 ± 0.019 


0.077 ± 0.020 



Table 4: Summary of results for the first eight moments. The first group of lines (labelled a Q (N,q 2 ,m 2 ) = 
(x n ~ 1 ) p qcd) gives perturbative moments for c and b production at the T(45') and Z° energies. The second 
group (labelled "Experimental data" ) gives the measured moments when explicitly given by the experimental 
Collaborations. In this case, the ISR correction (when applied) has been taken from the right panel of Fig. 3. 
The third group (labelled L>Np(iV) = (x Ar_1 )Np) gives the moments of the non-perturbative fragmentation 
function that we extracted from the data. The last three groups report the value of the parameters of the 
KLP, BCFY and PSSZ parametrization extracted from several data sets. 
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Figure 22: Values of a/m for D* and B mesons as a function of N. 



8.2 Implications for heavy-flavour hadroproduction 

It is legitimate to ask what is the impact of the new, high precision T(4S) 
data on the calculation of D meson production spectra in hadronic colli- 
sions, especially in view of the discrepancy between T(AS) and Z° data. The 
question has not, however, a straightforward answer. If the discrepancy is 
related to a power suppressed effect in the e + e~ coefficient functions, one 
should then privilege the Z° data, where power effects are much reduced. 
It is worth noting, however, that if we instead use the T(4S) data, the im- 
pact on the hadronic cross sections is quite limited. This is clearly visible in 
Fig. 19, where it appears that for N around 5 the T(45) moments are higher 
than the Z° ones by roughly 20%. This value is directly proportional to the 
D* production cross section in hadron collisions at large py Therefore, in 
this 'worst case' scenario, having used the ALEPH data (the only accurate 
ones available at the time), might have lead Ref. [29] to underestimate the 
D* hadronic cross section by 20%, an uncertainty which is anyway not larger 
than those of purely perturbative origin (variation of renormalization and 
factorization scales) or stemming from the parton distribution functions. 
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Figure 23: Ratios between new evaluations of the da/dpx production cross 
section of D* mesons in pp collisions at the Tevatron Run II and the value 
originally published in [29], (da/dpT) re {- 



These considerations are put on a more quantitative footing in Fig. 23, 
where we plot the ratios between new determinations of the px distribution 
of D* production at the Tevatron Run II and the central value obtained in 
Ref. [29]. The solid line, labeled 'BCFY, r = 0.1' is obtained by employing 
the same non-perturbative fragmentation function and the same parameter 
r = 0.1 as in [29]. Its small difference from one is essentially of perturbative 
origin. It is due to the different treatment of the perturbative fragmentation 
function in the FONLL code for heavy quark hadronic production [42,43] 
which, for consistency with the extraction of the non-perturbative parame- 
ters, has been modified to include also the Sudakov resummation for the ini- 
tial condition and the large- N regularization procedure described in Eq. (45). 
The five other curves are instead obtained with different non-perturbative 
forms and/or parameters relative to the T(4S) or to the Z° results from 
Table 4. As expected, using a different functional form (KLP) but a param- 
eter also extracted from the ALEPH data (a = 9.3) returns a result very 
similar to that of Ref. [29]. On the other hand, using determinations from 
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CLEO/BELLE data (a = 13) returns a larger cross section, the increase 
being of the order predicted above and not larger than the uncertainties of 
perturbative origin. 




N 



Figure 24: Non-perturbative moments from weakly decaying 5's data. 

As far as B mesons are concerned, the values for as are translated into 
non-perturbative moments in Fig. 24. The dotted band shows the values 
given by two extreme choices of olb- One can see that, using everywhere the 
value determined at N = 2, i.e. a# — 29, as done in [36] 10 , only overestimates 
the moment at N = 4 by a few percent. Up to iV = 8 the difference is never 
larger than 10%. Such an uncertainty is fully acceptable when calculating 
the hadronic production of B mesons, given the similar or larger size of 
the perturbative uncertainties and of those due to the parton distribution 
functions. 

10 In this reference, a pure NLL collinear resummation was used, without Sudakov re- 
summation and large-iV correction factor. This does not affect, of course, the small-N 
region, hence the determination of a very similar value for as- 
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9 Conclusions 



In the present paper, we have obtained two main results. First, we have 
shown that it is possible to perform excellent fits of D and B meson frag- 
mentation spectra in perturbative QCD, using all known results on the per- 
turbative heavy-quark fragmentation function, and compounding them with 
a simple parametrization of non-perturbative effects. For reasons of space 
we did not perform fits to available data on D s and A c production. We can 
provide the corresponding results upon request. 

A second striking result is the evidence of large non-perturbative effects, 
visible in the relation between the D* fragmentation function at the T(4S f ) 
and Z° energies. It would be interesting to understand the power law of 
these contributions. Their magnitude would suggest al/E scaling law. The- 
oretical arguments based upon infrared renormalons would favour, instead, 
a 1/E 2 behaviour. Because of the lack of precise D production data in the 
intermediate region, it is difficult, at this point, to discriminate between the 
two possibilities. We point out, however, that, if these non-perturbative 
corrections involve the coefficient functions, they may be present also in 
light-hadron production, where data at intermediate energy are available. It 
is thus possible that fits to the light-hadron fragmentation functions from 
T(4S) up to Z° energies may clarify this issue. 

The parametrization of the non-perturbative component of the heavy- 
quark fragmentation function is also relevant for the calculation of heavy- 
quark hadroproduction cross sections. In the present work we provide various 
related results, that can be used for such calculations. 
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